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ABSTRACT 

Using survey data, we have re-evaluated the correlation of flat spectrum radio sources with EGRET 
sources in the Northern sky. An analysis incorporating the radio and X-ray properties and the 7-ray 
source localization is used to gauge the reliability of associations and to search for counterparts of 
previously unidentified EGRET sources. Above |6|=:10°, where the classification is complete, wc find 
that 70% of the Northern EGRET sources have counterparts similar to the bright EGRET blazars. 
For several of these we identify known blazar counterparts more likely than the earlier proposed 3EG 
association; for ~20 we have new identifications. Spectroscopic confirmation of these candidates is in 
progress and we have found flat spectrum radio quasars and BL Lac counterparts with redshifts as high 
as 4. We also find strong evidence for a set of 28 objects with no plausible counterpart like the known 
EGRET blazars. These thus represent either a new extragalactic population or a population of Galactic 
objects with a large scale height. The survey has been extended into the plane, where we find several 
new blazar candidates; the bulk of the sources are, however. Galactic. Looking ahead to the GLAST 
era, we predict that several of the present 3EG sources are composite and that higher resolution data 
will break these into multiple blazar IDs. 

Subject headings: AGN: blazars - surveys: radio - surveys: optical ~ Gamma Rays 



1. INTRODUCTION 

The EGRET telescope on the GGRO satellite has de- 
tected 271 sources in a survey of the Gamma-ray (~ 
100 MeV to 10 GeV) sky. Of these roughly a quarter have 
been identified as blazars and along the Galactic plane 
there are a half dozen objects confirmed as young pulsars 
through their pulsed 7-ray emission. Thus the bulk of the 
sources remain to be identified. There is a young Galac- 
tic population along the plane clearly correlated with high 
mass stars (Kaaret & Cottam 1996; Yadigaroglu & Ro- 
mani 1997). An intermediate latitude excess, especially 
in the direction of the Galactic bulge, suggests an older 
Galactic population whose nature is uncertain. Finally 
there remain a substantial number of high latitude sources 
with no AGN identification. 

The selection of blazar candidates has largely proceeded 
by correlation with existing radio surveys (Hartman et al. 
1999; Mattox et al. 2001). In contrast, the Galactic plane 
sources and individual intermediate latitude sources have 
been the subject of targeted multi-wavelength campaigns 
(Roberts, Romani & Kawai 2001; Halpern et al. 2001; 
Wallace et al. 2002). In this project we attempt to ob- 
tain a more complete census of plausible blazar counter- 
parts, sifting sources with extant radio survey data and 
then conducting a multi-wavelength follow-up. We are ob- 
taining spectroscopic confirmation of the candidate AGN 
with Hobby»Eberly Telescope (HET) Marcario LRS spec- 
troscopy. This survey has already discovered a number of 
new likely 7-ray blazars, including good candidates for the 
most distant persistent 7-ray sources known. 

1.1. Blazar Properties 

The 'blazar' label is somewhat heterogeneous, but in 
the context of the unified AGN model, these sources are 
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believed to be viewed close to the axis of a powerful rel- 
ativistic jet. As such they are compact fiat spectrum ra- 
dio sources, with apparent superluminal motion at VLBI 
scales. The optical counterparts exhibit significant po- 
larization and OVV (optically violently variable) behav- 
ior (Urry & Padovani 1995). Optical spectroscopy often 
shows large equivalent width emission lines (fiat spectrum 
radio quasars) while a significant fraction are continuum- 
dominated BL Lac-type objects. The broad-band spec- 
tral energy distribution (SED) is sometimes used to di- 
vide these into two classes with 'red' blazars showing a 
synchrotron peak in the IR-optical with a Compton peak 
in the 7-ray while 'blue' blazars have a synchrotron com- 
ponent extending into the X-ray with a Compton peak 
inferred to extend to the TeV range (Urry 1999, and ref- 
erences therein). There are surveys underway to substan- 
tially increase the number of known blazars (e.g. DXRBS; 
Landt et al. 2001) which will eventually help in system- 
atizing the broad-band properties of these sources. 

Of these properties bright, flat spectrum radio emission 
seems best correlated with 7-ray activity, but this may 
be primarily a selection effect of present counterpart lists, 
which focused on the brightest radio sources as the prin- 
cipal candidates. In an effort to make a less biased census 
of the counterparts, we have re-examined the correlation 
with X-ray and radio properties in selecting candidates. 



2. COUNTERPART SELECTION AND 'FIGURE OF MERIT' 
RANKING 

The existing 7-ray blazar lists were largely selected from 
the Green Bank 6cm and 21cm single dish surveys (Con- 
don, et al. 1991; White & Becker 1992). With relatively 
poor resolution, confusion and extended jet/host galaxy 
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emission were a serious impediment to selecting flat spec- 
trum cores as blazar counterparts. This was especially 
true at lower radio fluxes and toward the Galactic plane. 
Accordingly few sub-Jy counterparts have been identified. 
Real progress can now be made since interferometric sur- 
veys covering much of the Northern sky are available at 
21cm (NVSS; Condon, et al. 1988) and 3.5cm (CLASS; 
Myers, et al. 2002 - the source list was kindly supplied in 
advance of publication by Ian Browne). The CLASS sur- 
vey, targeting conipac;t gravitational lens candidates, pre- 
selected flat spectrum sources by comparing NVSS 21cm 
and Green Bank single dish 6cm fluxes, mapping sources 
with spectral index a < 0.5 {Si, oc i^~") with the VLA 
A array, to resolve structure at the 0.2" scale. Coverage 
was complete for Dec> 0° and |5| > 10°; there was also 
partial coverage beyond these limits. We are primarily in- 
terested here in the isolated unresolved cores that are not 
lens candidates. 

By comparing the 3.5cm and 21cm interferometric fluxes 
we have a greatly improved estimate of the spectral index 
of the compact cores, which are well identified by sub- 
arcsecond positional matches. We find that a number of 
sources previously selected as 'flat' have extended high fre- 
quency emission. Of course with such a variable popu- 
lation, non-simultaneous observations can have erroneous 
spectral indices, but our estimates are certainly much bet- 
ter than extrapolations from the confused GB6 fluxes. 

When CLASS coverage is not available, wc do nec;d to 
use GB6 fluxes. Following the CLASS prescription, we 
subdivide the 6cm flux between the coincident 21cm NVSS 
sources to make spectral index estimates. We then extrap- 
olate with the 4.85GHz/1.4GHz spectral index to 8.4GHz 
to estimate a flux in the CLASS range. In general, this ex- 
tension was needed in the Galactic plane where a plethora 
of 21cm sources make spectral indices unreliable. How- 
ever, in a number of cases wc are able to show that no 
target meets the CLASS survey criterion in the EGRET 
error box. This will be used to rule out typical blazar- 
type counterparts. All sources classifled as blazars ended 
up having a CLASS detection, except 3EG J2016-I-3657. 

We augment this radio selection by using the convenient, 
relatively uniform RASS bright and faint source catalogs. 
Counterpart candidates are evaluated for suitability of op- 
tical follow-up with USNO-Bl.O/POSS photometry. Our 
analysis will use the positions of these counterparts in the 
'Test Statistic' (TS) maps of the Third EGRET {mG) cat- 
alog, which plot the estimated likelihood of the observed 
7-rays as a function of point source location. 

2.1. Radio/X-ray properties 

High frequency radio emission is clearly a good discrim- 
inant. However to make the least biased selection of coun- 
terparts, we compute the over-density of sources near high 
latitude ( |6| > 20°) gamma-ray detections in bins of radio 
flux, spectral index and x-ray flux. To do this, we com- 
pare the number of sources detected within a given TS 
probability contour with many random realizations of the 
sky. The latter were obtained by shifting the EGRET TS 
maps on the sky in 2 degree increments, computing the 
random source TS values and correcting for variation in 
sky area from the survey cuts. The source counts within 
a given confidence contour as a function of, e.g., 8.4GHz 



flux, were then compared to the random realizations after 
normalization by the effective sky coverage. We define the 
excess fractional source density associated with the 3EG 
source as 

n= 77 (1) 

for each flux bin. The distribution for 8.4 GHz flux within 

the 95% confidence region is shown in figure 1. A simple 
least-squares power law fit to the binned data then defines 
our excess source density function for evaluation of indi- 
vidual candidates. This exercise shows significant excess 
at fluxes well below the ~Jy limit considered in earlier 
analyses. Other authors have also recently concluded that 
fainter radio sources should be considered as plausible can- 
didates (WaUace et al. 2002). 
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8.4 GHz CLASS Flux Spectral Index 

Fig. 1. — The excess fractional source density witiiin EGRET 
95% error contours as a function of; 8.4 GHz (CLASS) flux (left), 
and 1.4GHz/8.4GHz spectral index (right). 

We have similarly measured the excess of flat spectrum 
sources as a function of spectral index. Here the corre- 
lation is strongest for a < —0.6, becoming negligible for 
a > 0.53. Correlation with the RASS source catalog also 
shows an excess. However the maximum excess for the 
brightest sources is only 0.5 and about half of the strong 
radio counterpart candidates were not detected in the rel- 
atively shallow all-sky survey (some are present in deeper 
pointed observations). Clearly we do not wish X-ray non- 
detection to exclude a candidate, although X-ray bright 
sources are somewhat more likely to be correlated with 
the EGRET detections. We thus add 0.5 to the best flt 
X-ray n so that it runs from 0.5 to 1.0. In this way X-ray 
detection can increase the source's selection by up to a fac- 
tor of 2 (for >1PSPC cps), but X-ray non-detection is not 
seriously constraining. Thus the X-ray component of our 
analysis is not crucial - if it is excluded no sources are lost 
from our final source list, although their detailed ranking 
does change somewhat. The excess probability functions 
exhibited no strong dependence on the gamma-ray source 
localization. Therefore, the area within the 95% contours 
were taken to be representative of the gamma-ray sky for 
the derivation of the merit functions. 

2.2. Source Position Weighting 

The combination of 8.4GHz flux, spectral index and X- 
ray flux provides a useful tool for identifying EGRET -like 
blazars. Indeed, we are studying a set of objects selected 
by these cuts, without any 3EG coincidence, as possible 
GLAST-detectable blazars. Given the strong 7-ray vari- 
ability and the limited duty cycle in the flaring state, it 
would not be very surprising for relatively bright sources 
to have eluded detection during the limited EGRET expo- 
sure. However, to select individual EGRET counterparts, 
we will also use the radio source position in the likelihood 
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analysis. Previous blazar identification efforts have used 
the circular or elliptical fits to the TS maps as an approx- 
imation to the allowed sky region. However, many TS 
maps contain open or complex contours at the 95% level, 
so we employ instead the precise TS at the radio source 
position. 

The value at our counterpart candidate position {a, 5) is 
compared to the map maximum, giving ATS = TSmax — 
TS{a, 5) . ATS has been calibrated by Monte Carlo sim- 
ulations (Mattox et al. 1996) to give the statistical proba- 
bility L[ATS{a, 6)] that the source lies within the contour 
of constant ATS for 50%, 68% 95% and 99% confidence 
levels. We linearly interpolate between these values to ob- 
tain a smooth L[ATS{a,S)]-, sources at large ATS* have 
little probability of association with the 3EG source, even 
if very blazar-like. The statistical error estimate for ATS 
does not, however, account for position uncertainties as- 
sociated with nearby unmodeled 7-ray sources and with 
other systematic and instrumental biases. This is con- 
firmed by plotting the distribution of bright 3EG blazars 
(Mattox et al. 2001) as a function of ATS. Considerably 
fewer than half of the sources lie within the 50% contour, 
and several lie outside the 99% contour. We find that the 
additional spread can be well modeled by dividing the map 
value of ATS by 1.5 and then computing the nominal sta- 
tistical positional probability. The probabilities renormal- 
ized in this way are a good match to the observed source 
distribution. 

2.3. Combined Figure of merit 

We combine the fractional excesses for each of the coun- 
terpart candidate's properties with the estimate of the po- 
sitional probability to compute a total 'Figure of Merit' 
for the source as a radio/X-ray counterpart: 

FoM = nsAGHz X ria X n^-ray X L{a, 6) (2) 
Notice that this product of excess source fraction and posi- 
tional uncertainty is not a normalized probability of source 
identification. It is however an unbiased ranking of the 
counterpart likelihood that can be compared across the 
entire (high latitude) 3EG population. To avoid confu- 
sion, we multiply by 100 in the quoted FoM values. 

To evaluate the significance of this FoM statistic, we 
produced random realizations of the radio/X-ray sky by 
random, independent RA and DEC draws from the true 
CLASS position list, after excluding the actual positions 
of the ^ 70 highest-FoM coincident sources [which would 
otherwise imprint the observed correlation on the random 
realizations] . We compute the average distribution of FoM 
for 10^ realizations of such random skies, and compare 
with the true FoM distribution in Figure 2. The hashed 
region shows the ±la range for our estimate of the fraction 
of sources in a given FoM bin that are 'real', i.e. in excess 
of random coincidence. For FoM< 0.1 there is little ex- 
cess correlation. There is a puzzling deficit of true sources 
with FoM^ 2, but we have not been able to trace this 
to any one of the FoM source properties and so conclude 
that this is a statistical fluctuation. We select sources with 
FoM> 0.25 (i.e. < 20% false positive, even for the low- 
est FoM) as good counterparts; this is about a factor of 
two above the lowest bins showing significant correlation. 
To facilitate comparison with earlier ID lists, we divide 
these sources in half with FoM> 1 designated as 'likely' 



counterparts and 0.25 < FoM < 1.0 as 'plausible' coun- 
terparts (note however that our 'plausible' sources have a 
relatively high fraction of true associations). The line in 
figure 2 shows a simple linear fit to the source probabil- 
ity with a strong decrease below FoM=0.1. Integrating 
through the distribution of measured FoM, we find that of 
the 35 'likely' sources we expect less than 3 false positives 
(> 92% good IDs) and of the 32 'plausible' sources < 6 
may be false positives (> 82% good IDs). We can compare 
the FoM distributions of our IDs with those designated in 
the 3EG catalog and in Mattox et al. (2001). Many of the 
highest FoM sources are common to all lists. However, 
both the 3EG and Mattox lists claim high confidence for 
several sources that we only assign low probability. In ad- 
dition each list has a number of sources that are not flat 
spectrum or compact and do not meet our interferomet- 
ric selection criteria. These sources, 11 in the 3EG and 7 
in Mattox et al. (2001) are mostly lower confidence ('a' 
and 'plausible' designations in these catalogs). However, 
in addition to selecting previous IDs, we also find ^ 50% 
more high confidence sources than either, and virtually all 
of our 'plausible' sources are new. Noting that Mattox 
et al. (2001) retained 'plausible' sources with estimated 
likelihoods as low as a few percent, while our lowest FoM 
sources have a likelihood > 80%, we believe that our list is 
less biased, more complete and more reliable than earlier 
efforts. In particular, we can more easily identify multi- 
ple associations with a given EGRET source, rather than 
taking only the 'best' blazar ID. 




Log(FoM) 

Fig. 2. — Random (histogram) and true (points, with Poisson 
error bars) distributions of our FoM. The fractional excess (true 
ID fra<;tion) for each bin is shown by the (Poisson) error range of 
the shaded region (rigiit scale). Note the rapid fall-off of reliable 
identification below Log(FoM)Ri —0.9. 

2.4. Radio Counterpart IDs 

We find at least one blazar candidate counterpart for 66 
of the Northern 3EG sources. In many cases, the previ- 
ously claimed 3EG 'high confidence' and 'plausible' blazars 
are recovered with the largest FoM value for a given error 
ellipse. However, in a number of cases, our criteria select 
a different blazar as much more likely than the claimed 
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catalog association. In addition, we find some 2 dozen 
new radio associations in the Northern sky, whose FoM 
confidence are at least as good as that of the lowest of 
the previously identified 'high confidence' sources. In sev- 
eral cases there arc multiple radio blazar IDs in a single 
EGRET error ellipse. One rather complex example is 3EG 
J0118-I-0248. In the 3EG catalog this was possibly iden- 
tified with 0119-1-041; our analysis does select it as a (low 
FoM) possible association, but it is well outside the 99% 
error contour. In the Mattox et al. (2001) analysis, the 
3EG source was attributed to 3C37 {z = 0.672); this AGN 
however has a quite steep spectrum for the compact com- 
ponent and docs not meet our blazar identification crite- 
ria. Instead, our analysis selects two compact flat spec- 
trum sources as the most likely counterparts. The first is 
a relatively low luminosity radio galaxy core at z=0.047, 
the second is a newly discovered flat spectrum radio QSO 
at z=4.0. The 3EG error contours are quite elongated 
with the major axis spanning these two sources (Figure 
3), although 0119+041 is also a plausible source of this 
extension. A second example is 3EG J0808-I-5114, which 
hosts a flat spectrum radio quasar at z=1.14 and a BL 
Lac at z=0.13, again separated along the major axis of an 
extended 7-ray uncertainty region. Given that the 3EG 
survey is strongly flux limited, that the 7-ray blazar lu- 
minosity function is quite steep, and that the poor 7-ray 
resolution causes substantial source overlap and confusion, 
it is not surprising that in some c;ases a combination of 
fainter sources can push a location above the detection 
threshold. In several other cases support for multiple IDs 
comes from complex or elongated likelihood ('TS') maps. 
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Fig. 3. — TS maps of possible composite 3EG sources. Left 3EG 
J0118-I-0248. The 3EG identification 0119-1-041, the steep spectrum 
Mattox et al. (2001) counterpart 3C037 (diamond) and our two 
new blazar counterparts (along the uncertainty region major axis) 
are shown. Right: 3EG J0808-t-5114. Again two high confidence 
identifications lie along the major axis. 

2.5. Radio non-IDs 

Since the GB6/NVSS/RASS provide a uniform survey 
(abetted by CLASS follow-up), the absence of plausible 
IDs for some EGRET sources is significant. For exam- 
ple, we find no IDs for two high latitude sources believed 
to be likely pulsars (3EG J0010-h7300= CTA 1 and 3EG 
J1835-I-5918, see Reimer et al. 2001). Our technique 
allows us to select individual 3EG sources as non-blazar 
IDs. This is particularly useful at intermediate latitude, 
as it lets us separate the blazar fraction of the popu- 
lation, leaving a cleaner non-blazar sample. North of 
Dcc=0°, 28 previously unidentified sources are selected as 
not like EGRET blazars. There is some concentration to 
the Galactic plane, but much of this sample is at high lat- 
itude, suggesting either a new class of extragalactic 7-ray 



sources or an old high scale height Galactic populations of 
pulsars or similar objects. 

2.6. Summary Results 

We summarize our source classification in Table 1 and 

Figure 4. In summary, of the 116 Northern 3EG sources 
(excluding the Solar flare), 66 have at least one plausible 
blazar- like radio counterpart. In the Table our high con- 
fidence associations are listed in bold face and our lower 
confidence associations in plain text. Previously claimed 
AGN associations not supported by our analysis are given 
in italics in parentheses, with FoM values computed when 
possible. The Table also notes 10 sources associated with 
pulsars or known plerion/young SNR (indented). These 
are neutron star/pulsar/SNR associations from the 3EG 
catalog, with a few additions from new discoveries which 
are noted in the individual source comments. Of course 
many additional associations have been claimed for Galac- 
tic objects; we do not evaluate or exhaustively review these 
here, but we believe that beyond the few pulsed detections, 
the identifications rapidly become quite speculative. 

An additional 28 sources are unidentified but accord- 
ing to our analysis are definitively unlike the known 3EG 
blazars. The remaining sources are at low |6| < 5° lati- 
tude. Only faint fiat spectrum sources are allowed in the 
error boxes of 7 of these and, given the high density of 
Galactic sources, these are almost certainly chance coinci- 
dence. The final 5 show bright extended HII regions which 
could, in principal, hide a bright flat spectrum compact 
blazar, but again in practice we expect the sources to be 
Galactic. Further radio imaging study could better rule 
out AGN counterparts in these cases. 

Figure 4 summarizes the present 3EG source classiflca- 
tions in Galactic coordinates. Above the Dec=0 line, we 
use our new classifications from Table 1, with blazars as 
filled circles. In the South, the same symbols are used for 
the 3EG 'A' and 'a' blazar designations. Pulsar IDs and 
strong pulsar /plerion candidates are shown by filled and 
open stars. The definitive non-blazars arc shown by open 
circles. The improved completeness of the Northern clas- 
sification is evident. Clearly pushing south into the bulge 
population will be very interesting. 

Since AGN are expected to be isotropic, it is of inter- 
est to examine the distribution of the northern blazar IDs. 
Detailed comparison requires folding through the EGRET 
exposure and sensitivity maps for a model blazar lumi- 
nosity function; we defer this to a future communication. 
However since exposure is relatively uniform for |6| > 20°, 
we report simple number counts. From 30° < |&| < 90° we 
find 36 blazar identifications, or 11.6/sr within our sur- 
vey boundary. Extrapolating toward the plane, we find 
for 20° < |6| < 30° 12 sources against 11 expected. At 
10° < |6| < 20° we find 8 sources vs. 12 expected, and a 
similar fraction within 10° of the plane. Given the higher 
background and detection threshold near the plane, this 
does not however mean that our identification of blazars 
at low |6| is incomplete. 
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3. OPTICAL FOLLOW-UP 

Of our radio-selected counterparts, 50 are new sources, 
not selected in previous EGRET coimtcrpart lists. Fur- 
ther, several additional selected radio sources, while 
flagged as possible counterparts in earlier studies, have 
not had spectroscopic confirmation of their AGN nature. 
We first checked for published spectroscopy of these radio 
sources, comparing with the Tenth Edition of the Quasar 
Catalog (Veron-Cetty et al. 2001) and cross checking by 
querying NED. We find that 16 (5 high probability IDs) 
of our newly selected counterparts have known redshifts. 
However over thirty sources in the Northern sky were 
found to have previously unstudied blazar-like radio/X- 
ray candidate counterparts, many with quite high FoM. 
These are the target of our HET spectroscopy, along with 
a number of similar sources having no iJCJiJiJT association. 

Just before submission, we learned that Halpern, Era- 
cleous & Mattox (2003) have recently completed follow-up 
spectroscopy of a set of previously claimed radio counter- 
parts of 3EG (and 2EG, GEV) sources with lower 4.8GHz 
fluxes. Their list includes 8 sources also identified in our 
exercises. For three sources we have new spectroscopy in 
common; for two their new redshift agrees with that which 
we obtain (below). We were also able to obtain a redshift 
for the third source (J1826-I-0149). For three other sources 
Halpern, Eracleous & Mattox (2003) obtain no z, for two 
the z they measure was already found in the QSO catalog, 
and for the last (J1239+0443) they obtain a new z. This 
value is included in our table as the only entry not from 
cataloged data or our new spectroscopy. 

3.1. HET LRS Spectroscopy 

Spectroscopy was obtained using the Marcario Low Res- 
olution Spectrograph (LRS)(HiU et al. 1998) on the 9.2m 
HET (Ramsey et al. 1998). This novel telescope, of tilted 
optical Arecibo design is still under development, but with 
a range of magnitudes the blazar candidates were suitable 
targets for the early operations phase. These targets were 
observed in regular queue observations from 3/02 -2/03. 
We obtained 2 x 450s exposures for most targets; a few 
of the most interesting fainter targets were observed with 
2 X 900s exposure. Observations were made employing a 
300 1/mm grating and a 2" slit for a dispersion of AA per 
(binned) pixel and an effective resolution of 16A covering 
AA4200 — 10500. The new spectra are displayed in Figure 
5. 

Standard IRAF CCD reductions, optimal extraction 
and calibration were performed. Unfortunately, the vary- 
ing pupil during the HET tracks is not yet fully moni- 
tored and so small (few %) fiat-field variations were not 
corrected, as with untracked flat field frames the varying 
illumination worsened the (already substantial) fringe fea- 
tures and sky subtraction errors in the near IR. Redshifts 
were measured by cross correlation analysis with AGN and 
galactic templates, using the IRAF RVSAO package. A 
few identifications are based only on the Mgll A2798 line, 
but in each case the proposed z is the only plausible value, 
given the absence of other expected strong lines. Several of 
our objects are continuum dominated with very faint line 
features, for these the substantial light grasp of the HET 
was important. All of the observed targets were found to 
be AGN, although some spectra were taken under rather 



poor conditions. Some of these have relatively uncertain 
z estimates, these are noted by (:) in Table 2. 
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Redshift 

Fig. 6. — The redshift distributions of our Northern blazar sample 
(solid line histogram), the Northern identifications from the 3EG 
catalog (shaded histogram) and the Mattox et al. (2001) sample 
(bar histogram). 

As of this publication we have obtained 29 new spec- 
troscopic identifications and an additional 16 new associ- 
ations with previous spectroscopy. The new and archival 
redshifts for our new blazar IDs are given in Tables 1 and 
2. The later table also contains precise positions and sky 
survey flux estimates; as expected for blazars, the flux 
at our observation epoch was often substantially different. 
Most of the AGN were found to be flat spectrum radio 
quasars. We denote sources with emission line equiva- 
lent width < 5A as BL Lacs. Several sources have the 
narrow lines (FWHM< lOOOkm/s) and lower excitation 
states characteristic of radio galaxies, we denote these as 
'R' in Table 1. It should be noted that within the classical 
AGN paradigm, narrow-line radio galaxies would not be 
expected to show blazar activity. The redshift distribution 
of our identified sources is shown in Figure 6. For compar- 
ison, we also plot the redshift distributions of the sources 
selected by the 3EG and Mattox et al. (2001) criteria. The 
z distributions arc rather similar, largely displaying spec- 
troscopic ID selection effects, but it is clear that we have 
substantially higher sensitivity to high redshift and have 
nearly doubled the maximum z value. 

4. NOTES ON INDIVIDUAL SOURCES 

3EG J0118+0248 As noted above, the Mattox et al. 
(2001) identification with the bright steep spectrum source 
3C37 is not supported by our analysis and the 3EG asso- 
ciation with 0119-1-041, is apparently superseded (or at 
least augmented) by our discovery of two flat spectrum 
candidates along the error ellipse major axis. If these con- 
tribute similarly to the 7-ray flux, they would represent 
extrema of the 7-ray blazar population, with a factor of 
~ 10^ difference in luminosity. If the association is con- 
firmed, J0122-F0310 at z=4.0 would be the highest redshift 
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non-GRB 7-ray source known. This 3EG source will be 
particularly interesting for study in the GLAST era, since 
with higher angular resolution the contributions of the sev- 
eral posited associations can be disentangled and GLAST's 
good sensitivity above lOGcV should allow intcrgalactic 
absorption of J0122+0310 against a cosmic background of 
starlight to be detected at > 20GeV. 

3EG J0204+1458 In addition to the hkely association 
J0204+1514 this source may also incorporate flux from the 
flatter and more centrally placed radio source J0205+1444. 

3EG J0215+1123 with a reasonably large FoM 
J0213-I-1213 is an interesting low z BL Lac candidate for 
TeV studies. 

3EG J0222+4253 This object has been argued to be 
a composite by Kuipcr at al. (2000) with a possible 
pulsed detection of the millisecond pulsar PSR J0218+42 
accounting for part of the sub-GeV flux. The remainder 
was attributed (also in the 3EG catalog) to the BL Lac 
3C 66A. This source, with ai.4/8.4 — 0.62, lies just be- 
low the spectral cut-off for inclusion in our FoM. How- 
ever, 2.3 GHz VLBA data are available (Fey & Chariot 
2000) which give a much flatter index for the compact 
component, which we adopt in Table 1. This does make 
3C 66A highly likely. CLASS analysis also selects the 
215mJy source J0223-I-4259 (3C 66B, a radio galaxy with 
a core+jet a,t z = 0.02) as a plausible counterpart. Con- 
fusion in the NVSS image prevents detection of the core, 
so using the 15GHz core flux measured by Jackson, et al. 
(1993) we obtain 08.4/15 = 0.38 and a strong plausible as- 
sociation. The GeV localization of Kuiper et al. (2000) 
supports either 3C object as the source of the higher en- 
ergy 7-rays. Certainly GLAST observations will be help- 
ful sorting out this complex region. Interestingly, a TeV 
detection is claimed in this region Neshpor et al. (1998) 
which would be somewhat surprising from z as large as 
0.44. The apparent position is consistent with cither 3C 
66A or 3C 66B and future air Cerenkov observations may 
be able to select the preferred source, given their 0.1° 
separation. 

3EG J0245+1758 The source closest to the TS maxi- 
mum is a strong- lined FSRQ at z=3.59. This source shows 
a likely damped Lya system at z « 3.15, raising the pos- 
sibility that the source is gravitationally magnifled by an 
intervening Galaxy. 

3EG J0329-I-2149 Mattox et al. (2001) noted the ra- 
dio source as a possible association. We confirm it as a 
likely blazar counterpart and have measured z=2.07, in 
agreement with Halpern, Eracleous & Mattox (2003). 

3EG J0404-I-0700 The better of the two likely associ- 
ations is our most puzzling case, spectroscopically. When 
first observed, the spectrum was highly continuum dom- 
inated. The two strongest lines were picked up as Lya 
and CIV A1449 at 2 = 3.13 in cross-correlation analysis. 
Both are (local) > Sc detections; they have rest equiva- 
lent widths of O.44A and 0.39A, respectively. CHI A1909 
is obscured due to the fringing and no other strong lines 
are expected. This redshift ID is supported by an appar- 
ent strong Mgll absorption doublet at z = 1.578 (REW 
O.I9/O.22I). With such low EW for the emission lines this 
source is (in this state) clearly a BL Lac (Marcha et al. 
1996). However a second observation several months later 
caught the source in a low state for which broad Mgll 



line and several Nc and O lines give a highly significant 
measurement of z = 1.133. The Mgll line is in fact also 
present in the high state at identical flux. The difference 
spectrum shows no Mgll line, but does have the higher z 
features noted above. We conservatively adopt the lower 
redshift, as z = 3.13 more than doubles that of the high- 
est z BL Lac cataloged to date (PKS 1309-216 z = 1.49 
Blades, Murdoch & Hunstead 1980). Further spectroscopy 
in a variety of flux states is certainly warranted. 

3EG J0450-f-1105 Halpern, Eracleous & Mattox 
(2003) have observed J0449-M121, identified in the 3EG 
catalog as the likely counterpart, and have not confirmed 
the catalog z. Their nearly featureless spectrum suggests 
a BL Lac. 

3EG J0459-I-3352 This Galactic error region has at 
least two maxima. The blazar candidate J0503-I-3404 is 
well associated with one peak, but the second peak, nearly 
2 degrees away is more likely Galactic. 

3EG J0634-I-0521 This source has been plausibly asso- 
ciated with the Be X-ray binary SAX J0635-h0533 (Kaaret 
2001), so it might be classified as 'p'. 

3EG J08084-5114 This source is probably composite 
with two likely counterparts (Figure 3). 

3EG J0917-I-4427 The bright, 1.3Jy, flat spectrum as- 
sociation is only plausible, being just outside the 99% con- 
tour. However, this error region is highly elongated and 
source is offset in the short direction. If an unmodeled 
systematic of ^ 0.3° induces this offset, this would be a 
likely association. 

3EG J1133-I-0033 and 3EG J1329-|-1708 we have 
observed a plausible radio source in each, identifying them 
as featureless BL Lacs, but have not obtained redshifts. 

3EG J1605-I-1553 Our new identification at z=0.11 
is brighter at 8.4GHz and has a smaller ATS than the 
association proposed by Mattox et al. (2001). 

3EG J1621-I-8203 Mukherjee et al. (2002) have previ- 
ously argued for identification with this Seyfert. 

3EG J1835-I-5918 This source has been the subject of 
intense study (Reimer et al. 2001; Halpern et al. 2002), 
who have both argued that it is an isolated Geminga-like 
pulsar. 

3EG J2016-I-3657 This crowded Galactic error region 
contains the SNR CTB 87, several bright HII regions in 
addition to the radio source B2013-I-370, claimed as a 
possible BL Lac (Halpern et al. 2001). Our spectrum of 
this source, taken under very poor conditions, is sufficient 
to support the continuum-dominated nature of the heav- 
ily absorbed counterpart, but is insufficient to determine 
a redshift. We tentatively adopt the blazar designation 
of these authors, but confirmation requires improved an- 
gular resolution and possibly variability correlation with 
GLAST. 

3EG J2021-I-3651 Is likely associated with the newly 

discovered pulsar and plerion discussed in Roberts et al., 
(2002); Roberts, Romani & Kawai (2001). 

3EG J2035-I-4441 May be associated with an X- 
ray/radio plerion (Roberts, Romani & Kawai 2001). 

3EG J2036-F1132 Our analysis produces FoM=0.22 
for B2032+107 (a z = 0.6 BL Lac), the 3EG A des- 
ignated counterpart of this source. This is just below 
our cutoff, but as it is significantly less than the FoM 
of J2034-I-1154, we propose re-assignment to this source. 
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Our second source J203H-1219 is bright and flat, but lies 
significantly outside the 99% confidence contour. A re- 
evaluation of the likelihood contours in terms of multiple 
sources could affect the rankings of these plausible coun- 
terparts. 

3EG J2206+6602 The 3EG catalog also selected this 
source with low probability; we confirm it as a z = 1.12 
FSRQ. 

3EG J2209+2401 This source was claimed as an iden- 
tification in (Mattox et al. 2001), but the listed rcdshift 
was evidently a typographical error, repeating that of an- 
other source. We have measured a single line z that ex- 
cludes the published result. 

3EG J2227+6122 This source is likely identified 
with the recently discovered X-ray/radio pulsar PSR 
J2229-I-6114 and its associated wind nebula (Halpern et al. 
2001b). 

5. CONCLUSIONS AND IMPLICATIONS 

We believe that we have made a relatively unbiased as- 
sessment of the association of flat spectrum radio sources 
with EGRET AGN, identifying plausible counterparts 
down to 8.4GHz fluxes of ~ lOOmJy. Assuming, as we 
have argued, that the bulk of our likely and plausible IDs 
are correct, we have substantially increased the complete- 
ness of identification of the Northern \b\ > 10° EGRET 
sources from ~40% to ~ 70%. We have also argued that 
in a number of instances the 3EG sources are composites; 
including this, the increase in the number of proposed 7- 
ray blazars is even larger. 

Our identification allows selection of sources with sub- 
stantially smaller 8.4GHz radio flux. It should not be 
too surprising that the number of candidates tapers off 
smoothly below the previous typical IJy value. We sus- 
pect that with smaller GLAST error ellipses or improved 
multi- wavelength constraints, identification could continue 
well below the ~ O.lJy limit of our survey. Our identi- 
fication of ~ 30 sources that have flat spectrum counter- 
parts (if any) well below O.lJy docs not, of course, preclude 
that some of these may be radio-faint or steep spectrum 
AGN. Indeed, excess coincidences do continue slightly be- 
low our FoM identification limit (Figure 2). Also a num- 
ber of the 3EG and Mattox et al. (2001) identifications 
not duplicated here are steep spectrum, but very bright 
low z AGN. In these cases (with the exception of 3EG 
J0118+0248/3C037, discussed above) the steep spectrum 
proposed counterpart is located at large ATS. Neverthe- 
less, if such sources are truly associated, it is reasonable 
that these form a distinct subset of the 3EG population, 
not identified in this analysis. However, the fact that an 
appreciable fraction of the 'non-blazar' 3EG sources cor- 
relate with the Galactic disk suggests that many represent 
a new (old) Galactic population. Detailed assessment of 
the completeness of blazar IDs through the plane requires 
careful treatment of the 3EG exposure maps; we defer this 
and other population analysis to a later paper. 

Along with fainter radio associations, we have pushed 



back the horizon of plausible blazar identifications (Figure 
4). It will be interesting to compare the SED and VLBI 
properties of these fainter blazar counterpart candidates. 
This survey has found several individual targets of partic- 
ular interest, including good 7-ray source associations at 
2 = 3 — 4, several distant BL Lacs, and several low z AGN 
that may be suitable targets for ground-based TeV ob- 
servations. These targets should be the blazars brightest 
to GLAST and they will provide the most detailed spec- 
tra and light curves; as such they merit further study in 
preparation for the GLAST era. Since spectroscopic iden- 
tifications are continuing, we defer detailed discussion of 
the red shift/ luminosity function distributions to a future 
paper. 

The 'non-blazars' are also excellent targets for further 
study in preparation for GLAST, since this sample will 
most likely produce new classes of high energy emitters. 
While we do not discuss here the correlation of 7-ray spec- 
tra and variability with the lower energy SEDs, we do 
note that our 'non-blazars' correlate fairly well with the 
'steady' sources of Gehrels, et al. (2000) and the 'persis- 
tent' sources of Grenier (2000), although exceptions oc- 
cur. Clearly, pushing the improved identifications south 
to cover the Galactic bulge and beyond will be very im- 
portant for characterization of these populations. 
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Fig. 4. — Aitoff equal area projection of 3EG sources in Galactic coordinates, showing our new classifications in the Northern hemisphere. 
Large filled circle=high confidence blazar, Smaller filled circle=plausible blazar, Filled star=pulsar, Open star=pulsar/plerion candidate, Open 
circle= Non-Blazar, cross=presently unclassified. Symbols south of DEC=0° are similar, with AGN drawn from the 3EG A/a classifications. 
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Fig. 5. — (continued) HET/Marccirio LRS spectra. 
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Note. — Column 7: Our new associations and/or redshifts * = new spectral ID. t= redshift from NED, 
10th QSO Catalog, etc. 

Column 8: Mattox et al. (2001) selected blazars +— 'High probability', — — 'Plausible' 
Column 9: 3rd EGRET Catalog blazars, A— 'High confidence', a—' lower confidence' 

Column 10: classification f— FSRQ, b— BL Lac, R— narrow-line radio galaxy, G— Likely galactic, N— 'Non- 
blazar', p=Pulsar candidate/plerion, P=confirmed pulsar. 
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Table 2 

New 3EG Counterpart Candidates with Spectroscopic Identification 



ID FoM a 5 R2 B2 z Type 

(J2000 Coordinates) Magnitude 
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05 
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+48 
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0.63 


05 
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08 
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+52 
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14.54 
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+00 


15 
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18.88 


19.43 


1.633* 




J1221+4411 
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Note. — Some sources have been previously flagged (see Table 1), but are spectroseopically 
confirmed here. Column: 5,6, USNO Bl.O magnitudes (Monet, et al. 2003), Column 7: new(*) 
or archival z, Column 8: classification (see Table 1). 



